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We examine metastable and transient effects both above and below the first-order disorder driven 
decoupling line in a 3D simulation of magnetically interacting pancake vortices. We observe pro- 
nounced transient and history effects as well as supercooling and superheating between the ordered 
and disordered phases. In the disordered supercooled state as a function of DC driving, reorder- 
ing occurs through the formation of growing moving channels of the ordered phase. We find that 
hysteresis in V(I) is strongly dependent on the proximity to the decoupling transition line. 
PACS numbers: 74.60Ge, 74.60Jg 



Vortices in superconductors represent an ideal system 
in which to study the effect of quenched disorder on elas- 
tic media. The competition between the flux-line inter- 
actions, which order the vortex lattice, and the defects in 
the sample, which disorder the vortex lattice, produce a 
remarkable variety of collective behavior [jlj . One promi- 
nent example is the peak effect in low temperature su- 
perconductors, which appears near H C 2 when a transition 
from an ordered to a strongly pinned disordered state oc- 
curs in the vortex lattice @-||. In high temperature su- 
perconductors, particularly BSCCO samples, a striking 
"second peak" phenomenon is observed in which a dra- 
matic increase in the critical current occurs for increasing 
fields. It has been proposed that this is an order-disorder 
or 3D to 2D transition. 

Recently there has been renewed interest in transient 
effects, which have been observed in voltage response 
versus time curves in low temperature superconductors 
In these experiments the voltage response in- 
creases or decays with time, depending on how the vor- 
tex lattice was prepared. The existence of transient states 
suggests that the disordered phase can be supercooled into 
the ordered region, producing an increasing voltage re- 
sponse, whereas the ordered phase may be superheated 
into the disordered region, giving a decaying response. In 
addition to transient effects, pronounced memory effects 
and hysteretic V(I) curves have been observed near the 
peak effect in low temperature materials. [7j. Xiao 
et al. Jt| have shown that transient behavior can lead to 
a strong dependence of the critical current on the cur- 
rent ramp rate. Recent neutron scattering experiments 
in conjunction with ac shaking have provided more di- 
rect evidence of supercooling and superheating near the 
peak effect jL5|. Experiments on BSCCO have revealed 
that the high field disordered state can be supercooled to 
fields well below the second peak line JrI. Furthermore, 
transport experiments in BSCCO have shown metasta- 
bility in the zero-field-cooled state near the second peak 



as well as hysteretic V(I) curves |17j and transient ef- 
fects p!s[ . Hysteretic and memory effects have also been 
observed near the second peak in YBCO [fi9|-pl|. 

The presence of metastable states and superheat- 
ing/supercooling effects strongly suggests that the order- 
disorder transitions in these different materials are first 
order in nature. The many similarities also point to a 
universal behavior between the peak effect of low tem- 
perature superconductors and the peak effect and second 
peak effect of high temperature superconductors. 

A key question in all these systems is the nature of 
the microscopic dynamics of the vortices in the transient 
states; particularly, whether plasticity or the opening of 
flowing channels are involved [^| . The recent experiments 
have made it clear that a proper characterization of the 
static and dynamic phase diagrams must take into ac- 
count these metastable states, and therefore an under- 
standing of these effects at a microscopic level is cru- 
cial. Despite the growing amount of experimental work 
on metastability and transient effects in vortex matter, 
these effects have not yet been investigated numerically. 

In this work we present the first numerical study of 
metastability and transient effects in vortex matter near 
a disorder driven transition. We demonstrate that the 
simulations reproduce many experimental observations, 
including superheating and supercooling effects, and then 
link these to the underlying microscopic vortex behav- 
ior. We consider magnetically interacting pancake vor- 
tices driven through quenched point disorder. As a func- 
tion of interlayer coupling or applied field the model 
exhibits a sharp 3D (ordered phase) to 2D (disordered 
phase) disorder-driven transition |22|. By supercooling 
or superheating the ordered and disordered phases, we 
find increasing or decreasing transient voltage response 
curves, depending on the amplitude of the drive pulse 
and the proximity to the disordering transition. In the 
supercooled transient states a growing ordered channel 
of flowing vortices forms. No channels form in the su- 
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perheated region but instead the ordered state is homo- 
geneously destroyed. We observe memory effects when 
a sequence of pulses is applied, as well as ramp rate de- 
pendence and hysteresis in the V(I) curves. The critical 
current we obtain depends on how the system is prepared. 

We consider a 3D layered superconductor containing 
an equal number of pancake vortices in each layer, in- 
teracting magnetically. We neglect the Josephson cou- 
pling which is a reasonable approximation for highly 
anisotropic materials. The overdamped equation of mo- 
tion for vortex i at T — is = — X^'i ^i^(Piji z ij) + 
i" p + id = Vj. The total number of pancakes is N v , and 
Pij and are the distance between vortex i and vortex j 
in cylindrical coordinates. We impose periodic boundary 
conditions in the x and y directions and open bound- 
aries in the z direction. The magnetic interaction energy 
between pancakes is |23|,|24[] 
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eo = *67(. 47r £.T> d = 0.005A is the interlayer spacing, A 
is the London penetration depth and £ = 0.01A is the co- 
herence length. When the magnetic field H increases, the 
distance p between pancakes in the same plane decreases, 
but the distance d between planes is unchanged. Thus we 
model H by scaling the strength of the in- and inter-plane 
interactions via the prefactor s m , such that s m oc 1/H. 
We denote the coupling strength at which the sharp 3D- 
2D transition occurs as s^,. We model the pinning as 
N p short range attractive parabolic traps that are ran- 
domly distributed in each layer. The pinning interaction 



is f- P = E£i(/p/&)(* " ^ } )e(fe - Jr, - i{ p) |)/A), 
where the pin radius £ p = 0.2A, the pinning force is 
f p = 0.02/q , and /q = e /A. Throughout this work 
we will use 16 layers in a 16A x 16A system with a vortex 
density of n v — 0.35/A 2 and a pin density of n p = 1.0/A 2 
in each of the layers. There are 80 vortices per layer, 
giving a total of 1280 pancake vortices. 

For sufficiently strong disorder, the vortices in this 
model show a sharp 3D-2D decoupling transition as a 
function of s m or H 22 . A dynamic 2D-3D transition 
can also occur |2^] . In the inset of Fig. |l|(b) we show the 
critical current f c and z-axis correlation C' z as a func- 
tion of interlayer coupling s m , illustrating that a sharp 
transition from ordered 3D flux lines to disordered, de- 
coupled 2D pancakes occurs at = 1.2. Here f c is 
obtained by summing V x = (1/N V ) v x and identi- 
fying the drive fd at which V x > 0.0005, while C z = 
1 - {{\r hL - r i)L+ i|/(ao/2)) 6(a /2 - \r hL - r ilL+1 \), 
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FIG. 1. Transient voltage response curves V x versus time, 
(a) Increasing response for supercooled state, for s m — 2.0 
and f d /fo = (bottom) 0.0053, 0.0055, 0.006, 0.0065, 0.007, 
and 0.0075 (top), (b) Decreasing response for superheated 
state, for s m = 0.7 and f d /fi = (bottom) 0.006, 0.007, 0.008, 
0.009, 0.010, and 0.011 (top). Inset to (b): Critical current 
f c (filled squares) and interlayer correlation C z (open circles) 
for varying s m in a system with 16 layers, showing the sharp 
transition from coupled behavior at s m < 1.2 to decoupled 
behavior at s m > 1.2. (c) Transient interlayer correlation 
response C z , for: s m = 2.0 and fd = O.OO6/0 (heavy line); 
s m = 0.7 and f d = 0.010/,* (light line). 

where do is the vortex lattice constant. The ordered 
phase has a much lower critical current, f° = 0.0008/q 
than the disordered phase, / c do = 0.0105/ *. 

To observe transient effects, we supercool the lattice by 
annealing the system at s m < s c m into a disordered, de- 
coupled configuration. Starting from this state, at t = 
we set s m > s c m such that the pancakes would be or- 
dered and coupled at equilibrium, apply a fixed drive fd 
and observe the time-dependent voltage response V x . In 
Fig. [l](a) we show V x for several different drives fd for a 
sample with s m = 2.0 in a state prepared at s m = 0.5. 
For fd < 0.0053/g the system remains pinned in a de- 
coupled disordered state. For fd > 0.0053/g a time de- 
pendent increasing response occurs. V x does not rise in- 
stantly but only after a specific waiting time t w . The 
rate of increase in V x grows as the amplitude of the fd 
increases. As shown in Fig. [l](c), C z exhibits the same 
behavior as V x . 

In Fig. 0(b) we show a superheated system prepared at 
s m = 2.0 in the ordered state, and set to s m = 0.7 at t = 
0. Here we find a large initial V x response that decays. 
With larger fd the decay takes an increasingly long time. 
The time scale for the decay is much shorter than the 
time scale for the increasing response in Fig. [l](a). In 
the inset of Fig. ^(b) we demonstrate the presence of a 
memory effect by abruptly shutting off fd- The vortex 
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FIG. 2. Pancake vortex positions (circles) and trajecto- 
ries (lines) for the supercooled phase at s m — 2.0 and 
fd = 0.007/o, seen from the top of the sample, at times: 
(a) t = 2500; (b) t = 7500; (c) t = 12500; (d) t = 20000. Pan- 
cakes on a given layer are represented by circles with a fixed 
radius; the radii increase from the top layer to the bottom. 
An ordered channel forms in the sample and grows outward. 

motion stops and when fd is rc- applied V x resumes at 
the same point. We find such memory on both the in- 
creasing and decreasing response curves. The response 
curves and memory effect seen here are very similar to 
those observed in experiments Q . 

In Fig. U we show the vortex positions and trajectories 
in the supercooled sample with s m — 2.0 from Fig. 0(a) 
for fd — 0.007/q f° r different times. In Fig. ||](a) at 
t = 2500 the initial state is disordered. In Fig. |](b) 
at t = 7500 significant vortex motion occurs through the 
nucleation of a single channel of moving vortices, which 
forms during the waiting time t w . Vortices outside the 
channel remain pinned. In Fig. ||(c) at t = 12500 the 
channel is wider, and vortices inside the channel are or- 
dered and have recoupled. The pinned vortices remain in 
the disordered state. During the transient motion there 
is a coexistence of ordered and disordered states. If the 
drive is shut off the ordered domain is pinned but remains 
ordered, and when the drive is re-applied the ordered do- 
main moves again. In Fig. |2|(d) for t = 20000 almost 
all of the vortices have reordered and the channel width 
is the size of the sample. Thus in the supercooled case 
we observe nucleation of a microscopic transport channel, 
followed by expansion of the channel. 

The vortex positions and trajectories for a superheated 
sample with s m = 0.7 and f d = 0.006/q , as in Fig. |l|(b), 
are shown in Fig. |](a-d). In Fig. |^(a) the initial vortex 
state is ordered. In Fig. y(b-d) the vortex lattice becomes 



50 o :lt r 
Q OiQ — q — © — O — © — © 

co — ® -4a© — & -%p — o — © — ^ 

© — © — © — © — — O — O .*. -a 
© — © . — O -^-O — O «<3 C-^h;, 

o — a — Q — o — <© -o -® ego 
o© © — o o — © o •• o 
o - a -e© —® — ® — © — ©: 

<© 0—3>^0 O—O 

O O ~i£> O Cfe^lg-© — © 

c@ ® oo o o 

(a) 



d2 ap 




(b) 



™9 o 



(c) " (d) 

FIG. 3. Vortex positions (circles) and trajectories (lines) 
for the superheated phase at s m — 0.7 and fd = 0.006/q , at 
times: (a) t = 500; (b) t = 1500; (c) t = 2500; (d) t = 4500. 
Each vortex line breaks apart into pancakes as it encounters 
the pinning sites. 

disordered and pinned in a homogeneous manner rather 
than through nucleation. Each vortex line is decoupled 
by the point pinning as it moves until the entire line 
dissociates and is pinned. 

We next consider the effect of changing the rate Sfd 
at which the driving force is increased on V(I) in both 
superheated and supercooled systems. Fig. [|(a) shows 
V x versus fd, which is analogous to a V(I) curve, for the 
supercooled system at s m = 2.0 prepared in a disordered 
state. V x remains low during a fast ramp, when the vor- 
tices in the strongly pinned disordered state cannot re- 
organize into the more ordered state. There is also con- 
siderable hysteresis since the vortices reorder at higher 
drives producing a higher value of V x during the ramp- 
down. For the slower ramp the vortices have time to 
reorganize into the weakly pinned ordered state, and re- 
main ordered, producing no hysteresis in V(I). 

In a superheated sample, the reverse behavior occurs. 
Fig. |(b) shows V(I) curves at different Sfd for a system 
with s m — 0.7 prepared in the ordered state. Here, the 
fast ramp has a higher value of V x corresponding to the 
ordered state while the slow ramp has a low value of V x . 
During a slow initial ramp in the superheated state the 
vortices gradually disorder through rearrangements but 
there is no net vortex flow through the sample. Such 
a phase was proposed by Xiao et al. [7| and seen in re- 
cent experiments on BSCCO samples [L6|. At the slower 
8fd, we find negative dVjdl characteristics which resem- 
ble those seen in low- flj^l and high- [^(| temperature 
superconductors. Here, V(I) initially increases as the 
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FIG. 4. (a): V(I) for a supercooled system at s m = 2.0. 
Light line: slow Sfd of O.OOOl/o every 2000 time steps. Heavy 
line: Sfd = 0.001/q. Inset to (a): Memory effect in the 
V x response of a superheated system with s m = 0.7 and 
fd = 0.009/q . The current is shut off for 400 time steps, 
(b): V(I) for a superheated system at s m ~ 1.1. From left to 
right, Sf d = (fast) 0.02/o* , 0.01/J, 0.005/ *, 0.002/ *, 0.001/ *, 
0.0002/o, and O.OOOl/o (slow). Inset to (b): Effect of super- 
cooling on f c . Filled squares: equilibrium f c . Open diamonds: 
f c for samples prepared in a disordered, decoupled state. 

vortices flow in the ordered state, but the vortices decou- 
ple as the lattice moves, increasing f c and dropping V(I) 
back to zero, resulting in an N-shaped characteristic. 

To demonstrate the effect of vortex lattice disorder on 
the critical current, in Fig. 4 we plot the equilibrium f c 
along with f c obtained for the supercooled system, in 
which each sample is prepared in a state with s m = 0.5, 
and then s m is raised to a new value above sf n before 
f c is measured. The disorder in the supercooled state 
produces a value of f c between the two extrema observed 
in the equilibrium state. Note that the sharp transition 
in f c associated with equilibrium systems is now smooth. 

Our simulation does not contain a surface barrier 
which can inject disorder at the edges. Such an effect 
is proposed to explain experiments in which AC current 
pulses induce an increasing response as the vortices re- 
order but DC pulses produce a decaying response |||||. 
We observe no difference between AC and DC drives. 

In low temperature superconductors, a rapid increase 
in z-direction vortex wandering occurs simultaneously 
with vortex disordering |L5]|, suggesting that the change 
in z-axis correlations may be crucial in these systems as 
well. Our results, along with recent experiments on lay- 
ered superconductors, suggest that the transient response 
seen in low temperature materials should also appear in 
layered materials. 

In summary we have investigated transient and 
metastable states near the 3D-2D transition by super- 



cooling or superheating the system. We find voltage- 
response curves and memory effects that are very similar 
to those observed in experiments, and we identify the 
microscopic vortex dynamics associated with these tran- 
sient features. In the supercooled case the vortex motion 
occurs through nuclcation of a channel of ordered mov- 
ing vortices followed by an increase in the channel width 
over time. In the superheated case the ordered phase 
homogeneously disorders. We also demonstrate that the 
measured critical current depends on the vortex lattice 
preparation and on the current ramp rate. 
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